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Abstract 

In this paper, we first apply the Fitzpatrick algorithm to osculatory rational 
interpolation. Then based on Fitzpatrick algorithm, we present a Neville-like 
algorithm for Cauchy interpolation. With this algorithm, we can determine 
the value of the interpolating function at a single point without computing 
the rational interpolating function. 
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1. Introduction 

Interpolation is an important method in numerical approximation. Ratio- 
nal functions sometimes are superior to polynomial for interpolating data be- 
cause they can achieve more accurate approximations with the same amount 
of computation In addition, rational interpolants have a natural way of 
interpolating poles whereas polynomial interpolants do not. So how to solve 
the problem of rational interpolation is what people have been concerning. 
There are rich literatures on the univariate Cauchy interpolation and 



osculatory rational interpolation problem, such as [2| 



mterpoiation ana 

. a Baa a a. 

For multivariate rational interpolation j2|, y, \% llO|, uJj gave some results 
about bivariate cases and the authors assume the interpolation nodes are 
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on rectangular girds. 12|] computed rational interpolation over pyramid- 
typed grids in by branched continued fractions. When the interpolation 
data are scattered, 13j presented a fast solver for the linear block Cauchy- 



Vandermonde system that translates the interpolation conditions. [IJ] used 
the theory of algebraic geometry to study the minimal multivariate rational 
interpolation. 

One of the main problems of rational interpolation is the parametrization 
of all solutions of a given degree of complexity. Based on Euclidean algorithm, 



15| investigated a general frame work which lead to a parametrization of all 



16| considered the set M = {(a, 6) 



rational interpolation functions. 
bh mod x^*} of all solutions of the key equation for alternant codes, and 



give the Fitzpatrick algorithm. [17| extended the Fitzpatrick algorithm to 



determining a parametrization of all minimal complexity rational functions 
a{x)/b{x) interpolating an arbitrary sequence of points, and complexity is 
measured in terms of max{deg(a(x)), deg(6(x)) + where ^ is a given 
integer, [isj presented an algorithm to seeking the Grobner basis for the 
solution of polynomial congruences in one or more variables, 
the work in 
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19[ generalized 



and got a general algorithm applicable to a wide range of 
constrained interpolation. 

In this paper we apply the Fitzpatrick algorithm which appears in 
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to multivariate osculatory rational interpolation, and get the parametric so- 
lution of the multivariate osculatory rational interpolation function r(x) = 
a{X)/b{X). Based on Fitzpatrick algorithm, we present a Neville-like algo- 
rithm for multivariate Cauchy interpolation. 

The rest of the paper is organized as follows. In Section 2, we present 
the Fitzpatrick algorithm. In Section 3, we apply the Fitzpatrick algorithm 
to seek the weak solution {a{X),b{X)) of multivariate osculatory rational 
interpolation. In Section 4, based on Fitzpatrick algorithm, we give a Neville- 
like algorithm for multivariate Cauchy interpolation. With this algorithm, 
we can determine the value of the interpolating function at a single point 
without computing the rational interpolating function. 



2. Fitzpatrick algorithm 

Fitzpatrick algorithm, also called FGLM-like algotithm, is applicable to 
coding theor y, F ade approximation, partial realization, interpolation, and 
other fields. |20| described these literatures from a historical point of view, 
started from [l6|, and covered recent developments for list decoding. For 
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further details, please refer to 20| or 2lj. Now we introduce the structure 
of Fitzpatrick algorithm. 

Let F be a field, V = ¥[ polynomial ring and > 1 be a 

natural number. 

We denote by {ei, . . . , e^} the canonical basis of V^. Any term in 
is of the form m = (pe^, 1 < k < d, where is a term in V, and the 
set of terms in is denoted by T^'^\ Let -< be a term order. For each 
/ = ^ c(/, r)r G V^, its support is 

supp{f):= {rGr'^:c(/,r)^0}, 

its leading term is LT(/) := max^{supp{f)) , its leading coef ficient is 
LC(/) := c(/,LT(/)) and leading monomial is LM(/) := LC(/)LT(/). 

Let Mfc, Mfc+i be submodules of a P-module M, with Mk 5 M^+i, such 
that, for each s, 1 < s < n, there exists /3s G F satisfying 

(x, - /3,)Mfc C Mfc+i. (1) 

For each k, there exists an F-homomorphism 

Ok-.Mk^W (2) 

with kerfi9fc) = Mk+i. 

In Q, they described Fitzpatrick algorithm in the following theorem. 



Theorem 1. 19|] Let M be a V -module and let Mk ^ Mk+i be submodules of 
M satisfying ([I]) and for suitable f3s, Ok- Let H : — )■ M be an ¥ -linear 
function such that for each s, 1 < s < there exists 7^ G F satisfying 

H{xsb) = {xs + is)H{b) 

for all b = (61, ... , bd) G V^. Let S C be a submodule satisfying 

H{b) = mod Mfc, ybeS (3) 

and let S' C S be the set of elements satisfying 

H{b) = mod Mk+i (4) 

Then S' is a submodule ofV^. 
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// we have obtained an ordered minimal Grobner basis G = {G[l],. . . , 
G[\G\]} of S with respect to a term order -<, then a Grobner basis G' of S' 
with respect to -< can be constructed as follows: 
Define a, = ek{H{G[j\)), for 1 < j < IGI- 
// aj = for all j then 

G' = G 
otherwise 

j* = the least j for which aj ^ 

Gi := {G[j] : J < f} 

G2 ■■= {{xs - i^s + ls))G[r] ■.l<s<n} 

Gs:={G[j]-{a,/a,,)G[f]:j>f} 

G":=GiUG2UG3. 

In the following section we will use the Fitzpatrick algorithm to compute 
osculatory rational interpolation functions. 

3. Osculatory rational interpolation and Fitzpatrick algorithm 

Now, we introduce some notations. Let a = (ai, ■ ■ ■ , a„) G N". We 
define a differential operator D"' by 

Qai-i hofn ^ 5"!^ 



ail - --aJ dxi^ ■ ■ ■ dx'^" a\ dx"^ ■ ■ ■ dx^"' 

A subset A C N" is called a lower set, if it is closed under the division order, 
that is, if a = (ai, . . . , a„) G A then /3 E A for all (3 = . . . , with 
/3i < ai,i = 1, . . . ,n. 

The multivariate osculatory rational interpolation problem can be stated 
as follows: 

Given a set of L distinct points {Yi, . . . , Yl} in space F". Point Yi has 
multiplicity defined by the lower set Ai, and the corresponding values {f- "^ G 
F : Va G Ai, i = 1, . . . , L}. Construct a rational interpolation function 

a(X) 

such that 

D'^r{X)\y = f^''^ , yaeA, t = l,...,L, 

where a e V, b e V, a 0, b{Yi) 7^ for all i. 

We call this problem multivariate osculatory rational interpolation, and 
r{X) multivariate osculatory rational interpolation function. 
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3.1. Weak interpolation 

From the definition of multivariate osculatory rational interpolation we 
know that the equivalent definition is: the Taylor series expansion of r{X) 
at the point X = Yi satisfies 

t{x) = 5^ (X - r.)"'-7^'^') + ■ ■ ■ . 

L 

Let Si = ^Ai, i = 1, . . . , L, N = ^Sj. For each point Yi and the lower set 

i=l 

Ai, define polynomial hi 

hr.= Yl ix-Yir^-fi^-^\ 

For each Ai, rearrange the elements of Ai, such that each subset Aij = 
{ai^o, . . ., aij}, < j < Sj — 1, is still a lower set. In particular, Ai,si-i = Ai. 

Denote the ideal I{{Yi,Aij)) = {p e V : Z^XX)]^ = 0,Va G Aij} by 
lij and call lij the vanishing ideal of (Yi,Aij). 

Definition 1. (Weak interpolation) A pair {a,b) G is called a weak 
interpolation for multivariate osculatory rational interpolation problem if 

a = bhi mod li^s^^i, i = 1, . . . , L. 

Define (a, b) + (c, d) = {a + c,b + d), d{a, b) = {da, db). Thus M = {{a,b) : 
a = bhi mod Ii,si~i, i = 1, . . . , L} is a. P-submodule. 

If {(oi, . . . , (flj, is a Grobner basis of the submodule M, then any 
pair (a, b) with the form 

(a, 6) = ci(ai,6i) H h Ct{at,bt) 

is a weak interpolation, where Cj E V [j = 1, . . . ,t) are free parameters. 
Choose Cj properly such that b{Yi) 7^ 0, i = 1, . . . , L, then we can get the 
interpolation function 

a{X) Ciai + ■ ■ ■ + 
b{X) ~ cibi + ■■■ + ctbt ■ 
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3.2. Fitzpatrick algorithm_RI 

The aim of this subsection is to apply the Fitzpatrick algorithm to com- 
pute osculatory rational interpolation. 

Definition 2. (order -<g) 

1. We say 0) X^{1, 0) if |a| < or |a| = and X" -<fex , 

2. We say 0) ^5 X^(0, 1) if |a| < + 

3. We say X°(0, 1) X'^(0, 1) if |a| < or |a| = \(3\ and X" ^z^x , 
where -<iex is the lexicographic order on "P, and ^ is a given integer. 

It is easy to check that the order is a monomial order on V^. 

For each Aij, l<i<L,0<j<Si — 1, define the congruent equation 

as 

a = bhi mod Jjj, 

where lij is the vanishing ideal of {Yi,Aij). 

Define an order < on the lower sets {Aij, 1 < i < L,0 < j < Si — 1} such 
that ^iiji < Ai2j2 if and only if ii < i-^ , 01 ii = 12 = i and Aij-^ C Ajjj 

3 1 < J2- 

Consequently, an order on the congruent equations is induced: 
a = bhi mod Jj^j-^ < a = bhi mod li^j^ if and only if Ai-^^j-^ < Ai^^^^. 
Now we can establish a one to one correspondence between index k and 
Define a sequence of submodules M^, k = 0, . . . , N, where Mq = V^, 
Mfc is the set of common solutions of the first k congruent equations 

Mfc = {(a, b) G Mfc„i : a = bhi^ mod li^j^, },k = 1, . . . ,N. 

Obviously M = M^q = {{a,b) : a = bhi mod Ii^si-i,i = 1, • • • ,L}, and 
Mo 3 Ml D ■ ■ ■ D Mtv. 

Fix an order we compute the minimal Grobner basis ^tv of re- 
cursively. It is easy to see that {(1,0), (0, 1)} is a Grobner basis of Mq. We 
compute the Grobner basis Qk+i of Mk+i through the minimal Grobner basis 
Qk of Mfc. 

Let Qk = {(ai, 61), . . . , (0^^., bm,.)} be the minimal Grobner basis of M^, 
and the {k + l)-th congruent equation be a = bhi mod I^kr Then 

1. if ki = 0, that is Ai^ki = Ai^ = {«/ = 0}, then 

bhi-a = v = v{X - y;)"''0 mod /;,o. 
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2. if ki 7^ 0, that is Ai^k\A,ki-i = then 

bhi-a = u'{X-Yir^^^i mod Ii^kr 

Therefore for any (a, b) G Mk, we have 

bh,-a = iy{X-Yi)'"-''i mod 

and (a, b) G M^+i if and only if = 0. 
Define an F-homonorphism 

9 : Mk — ^F 
(a, b) I — > V 

Obviously, ker(^) = M^+i, (x, - |/z,,)Mfc C M^+i. We define R: H{{a,b)) = 
{a,b). Then for any (a, 6) G V"^, H[{xs - yi,s){a,b)) = (x^ - yi^s)H{{a,b)). 
Let ^ = Mfc, y = Mfc+i. 

With the definitions above, we can give the form of Fitzpatrick algorithm 
for multivariate osculatory rational interpolation. 

Fitzpatrick algorithm_RI : (Using the minimal Grobner basis Qk, com- 
pute the minimal Grobner basis of M^+i) 

Input: the minimal Grobner basis Qk = {(«!, ^i), • • • , {cLm^, bmk)}', 
Output: the minimal Grobner basis of M^+i, that is Qk+i] 

Rearrange the elements of Qk such that LT(ai, bi) ■ ■ ■ LT(amfe, ^mfc); 
for t from 1 to rrik do 

bthi-at = i^tiX -Yir^^^ mod Ii,k,; 

end do; 

a Ut = for all t then 

Gk+i '■= Gk', 

else 

for t from 1 to do 

Find the least tk such that z/j^ 7^ 0; 
end do; 

for t from + 1 to do 

{at,bt) := (at,fet) {at^,btj; 
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end do; 

Gk+i ■= . . . , {at^-i,bt^_i),{at^,btj ■ {xi - yi^i), 

{at„btj ■ {xn-yi,n),iat^+i,bt^+i), . . .,(a^,,6„J j 

Qk+i '■= minimal Grobner basis(^fc+i) 
end if; 
return Qk+i, 

We must point out that here we do not require LC((a, b)) = 1 in minimal 
Grobner basis. 

Example 1. Given the interpolation problem 



point 




d £ 
dx J 


d_f 

dyJ 


ML f 

dxy J 


(-1,2) 


2 








(1,1) 


3 








(2,1) 


4 


5 


2 




(3,2) 


3 


4 


3 
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Tabic 1: interpolation 



Fix the order -<o, using the Fitzpatrick algorithm_RI, we can compute 
the minimal Grobner basis of the submodule M: 

(n h \ = ( -r^ ^^^'^ Til -J/^ I ^'^^^ T II 1 61 2 , 3047 , 

{Ui,uij 14528 y 7264 ^ 14528 7264 y ' 908 ^ 14528 



731 _ 6335 \ 
14528 y 14528/' 



nn h.) = ( 19899 2 I 43619 1999 2 153069 , 14059 , -i 122 , 

"2, "2; 314176 ^ 314176 y 157088^ 314176 ^ 157088 ^ ^ ' 4909 ^ ^ 

„,2 67507 ^ 19127 „, , 135787 \ 



19636^ 314176 314m ^^+gl4m _ 18057^ , , ^ JS^o , 

V"3,"3; V 371696 ~'~ 371696 ^ 185848 ^' 28592 ~'~ 185848 ^ ~'~ ' 23231 ~^ 
61 _ 89291 _ 12399 , 141603 \ 

1787 ^ 371696 371696 ^ ~'~ 371696 / ' 

(n. hA = - 305 ,,3 _ 1519 2 , 2959 , 1481 2 _ 10769 _ _673_ , 

Vu.4,<^4; V 12438^ 24876 ~'~ 24876 ^ ~'~ 12438^ 24876 6219^ ~^ 

1 i22_T7/ -^ -^7/2 _ 4847 _ 1697 , lOO^N 
' 6219 y ~'~ 1382|'g 24876^ 24871^ ~'~ 248J_(3^/ ' 

(0-5, ^5) = (988^2/ ~ 494^ ~ T^'^y ~^ 247^ 247*^ ~ V ~ 247' ^ 247'^^ 
988 y ~'~ 247 247 y 247/' 

(oe, be) = {§iX^y - jk^V^ + §fy^ ' §7^^ + ik^V - v'^ - %^ + + 

184 1 _ 5^ 2 _ 1^ _ , 115 Y 

237' 6 y 474^ 3 474^^ 237/' 

yUT,Uj) — y^^^X 326-^ y ^ 489'^y ^163-^ 978-^i' 163^ 978"^ ^ ^ ^ 
92_ ^^2__55_ 176 , _55_ , 115\\ 

163' 978 978 y 489 ~'~ 326^ ~ 326/ J 

Any weak interpolation (a, b) have the form 



(a, b) = ci(ai, 6i) H h 07(07, 67), 
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where Cj E V (j = 1, . . . , 7) are free parameters. Choose Cj properly such 
that h{Yi) 7^ 0, 2 = 1,...,//, then we can get the interpolation function 



a{X) ciOi + ■ ■ ■ + c-jaj 
h{X) ~ cihi + ■ ■ ■ + cjh' 

4. Neville-like interpolation 

Neville's algorithm is used for polynomial interpolation which was derived 
by Eric Harold Neville. The algorithm aims at determining the value of the 



interpolating polynomial at a single point x. [22| also derived a Neville type 
algorithm for univariable rational interpolation. 

In this section, we present a Neville-like algorithm for multivariate Cauchy 
interpolation based on Fitzpatrick algorithm_RI. 

Given a set of L distinct points {Yi, . . . , Y^}, Yj G = 1, . . . ,L, and 
the corresponding values {/i,...,/^}, fj G F, j = 1,...,L, we want to 
determine the interpolating value at the point lo- 
in this case, we know that hj := fj, and := {(a, b) G Mfc_i : bh^ — a = 
mod Ik} = {(a, 6) G Mk-i : (bhk - a)|^^ = 0}, A; = 1, . . . , L. In Fitz- 
patrick algorithm_RI, if we can get the values W{i,j) = {bihj — ai)\y = 
and y^y bi\^^^ y^^) recursively without computing the weak interpolation 

(a, b), then we can determine the interpolating value at the point Yq. It means 
that using the present values, we can calculate the new values of W{i,j) and 
^^*l(a;o yo)' ^*l(a;o 2/0)-^ without Computing the weak interpolation (a, fe) when a 
new point is added. Based on this idea we get a Neville-Like algorithm for 
Cauchy interpolation. For simplicity we will restrict ourselves to the case 
n = 2. Three and higher dimensional cases can be treated similarly. 

Fix the order -<^, and y x. Let {(«!, . . . , (0^^., bm,.)} be the minimal 
Grobner basis of M^. Define W{i,j) = {bihj — aj)|^, i = l,...mfc, j = 

1, . . . , L, W{^, L + l) = (a,| k\^^^^J, W{^, L + 2) = LT((a„ 6,)). 

Define W, = {W{t,l), . . . ,W{i, L),W{i, L + l),W{i, L + 2)). 
We know that ((1, 0), (0, 1)) is the Grobner basis of Mq = P^. 
First, using the Grobner basis of Mq, we compute 

W{l,j) = 0-hj-l = -l, j = l,...,L; 

W{2,j) = l-h,-0 = fj, j = l,...,L. 

Let 

W{1,L + 1) = (1,0), iy(l,L + 2) = LT((1,0)) = (1,0), 
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W(2, L + l) = (0, 1), W{2, L + 2) = LT((0, 1)) = (0, 1). 
We will compute the W^-*^'' recursively by Wi, W2. 
Define the matrix 

Let Wfc_i = (Wi, . . . , W^j,)-^. Denote by tlVVfc_i the number of rows. 

Neville-like algorithm: 

Input: Wo = (Wi,W2)^; L; 

Output: Wl; 

for k from 1 to L do 

rrik ■■= tl>Vfc-i; 

rearrange the elements of Wk-i so that 
Wfc_i = (Wi,...,W„J^ 
and W{1, L + 2) <^ WijUk, L + 2); 

for i from 1 to m^, do 

Find the least zq such that W{iQ, k) 7^ 0; 

end do; 

if W{i, k) = for all i then 

m ■■= m-i 

else 

for i from io + 1 to rrik do 
for j from 1 to L do 

W{t,j) := W{t,j) - W{to,jy, 
end do; 

W{i, L + l):= W{i, L + 1)- ^^j^ Wito, L + 1); 

W{i,L + 2) := Wit, L + 2); 
end do; 

for j from 1 to L do 

W{mk + 1, j) := W{io,3) ■ {xj - x^J; 
W{mk + 2, j) := W{io,j) ■ {yj - yj,); 
end do; 

W{mk + 1, L + 1) := W{io, L + 1) ■ {xq - x^J; 
W{mk + 1, L + 2) := W{io, L + 2) ■ x; 
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W>(mfc + 2, L + 1) := W{i^, L + 1) ■ (yo - l/io); 
H>(mfc + 2, L + 2) := (zq, + 2) ■ 
Wfc := (Wi, . . . , Wi„_i,W,,+i, . . . , W„,+2)^; 
Wfe := Minimal Grobner basis (Wfc) 
end if; 
end do; 
return Wl; 

Minimal Grobner basis(W) 
Input: W;_ 

Output: W = (Wi,...,W^)^ with no W{i,L + 2) is divisible by 
H>(j,L + 2)forz^j. 

Example 2. Given objective function ln(x^ + y^), we will use the values 
at the points (1.75,1.75), (2.25,1.75), (1.75,2.25), (2.25,2.25), (1.85,1.85), 
(2.15, 1.85), (1.85, 2.15), (2.15, 2.15) to estimate the value at the point (2, 2). 
Fix the order -<o, L = 8, Neville-like algorithm outputs Ws: 



0. 


0. 


0. 


0. 


0. 


0. 


0. 


0. 


(8.200 X 10-^ 4.100 X 10-^) 


(x^0) 


0. 


0. 


0. 


0. 


0. 


0. 


0. 


0. 


(12.45851102, 5.991319801) 




0. 


0. 


0. 


0. 


0. 


0. 


0. 


0. 


(-37.67095363, -18.11585643) 


(0,x?/) 
(0,x2) 


0. 


0. 


0. 


0. 


0. 


0. 


0. 


0. 


(37.67095324, 18.11585624) 


0. 


0. 


0. 


0. 


0. 


0. 


0. 


0. 


(-1.225410500, -0.5894116890) 


(l/',0) 


0. 


0. 


0. 


0. 


0. 


0. 


0. 


0. 


(3.470473588, 1.668910238) 


(V,0) 



Table 2: Wg for \n{x^ + y"^) 



From the (8 -|- l)-th column of Wg, we can see that each of the vectors 

-> II '^j I (2 2) 

W{i, 8 + 1) = (cj , , , 6j , ,), z = 1, . . . , 6, gives an approximate value — - — — 

1(2,2) l(2,2J Oil (2,2) 

of ln(22 + 22). 
Here we choose 













Yo 



as our estimation value(see Table |3]), where Yq = (2,2), sgn(a;) satisfies: if 
X > 0, sgn(x) = 1, else sgn(x) = -1. Actually ln(22 + 2^) = 2.079441542. 
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i 




\ii{xl + yl) 


interpolating value of Wj 


1 


(1.75,1.75) 


1.812378756 


1.312378756 


2 


(2.25,1.75) 


2.094945728 


1.812378756 


3 


(1.75,2.25) 


2.094945728 


2.122484930 


4 


(2.25,2.25) 


2.315007613 


2.107686660 


5 


(1.85,1.85) 


1.923518459 


2.082067864 


6 


(2.15,1.85) 


2.085050780 


2.082067864 


7 


(1.85,2.15) 


2.085050780 


2.079431546 


8 


(2.15,2.15) 


2.224082865 


2.079439873 



Table 3: estimation value of ln(22 + 2^) 



i 




v/1 - (x.)2 - (y,)2 


interpolating value of Wj 


1 


(0.45,0.45) 


0.7713624310 


0.6713624310 


2 


(0.55,0.45) 


0.7035623640 


0.7078673362 


3 


(0.45,0.55) 


0.7035623640 


0.7035623636 


4 


(0.55,0.55) 


0.6284902545 


0.7035623639 


5 


(0.5,0.45) 


0.7399324293 


0.7035623639 


6 


(0.5,0.55) 


0.6689544080 


0.7047928585 


7 


(0.45,0.5) 


0.7399324293 


0.7071486038 


8 


(0.55,0.5) 


0.6689544080 


0.7071187945 



Table 4: estimation value of ^1 - (0.5)2 _ (g 5)2 



Example 3. Given the values of yl — x^ — at the points (0.45,0.5), 
(0.55, 0.45), (0.45, 0.55), (0.55, 0.55), (0.5, 0.45), (0.5, 0. 55), (0.45,0.55), (0 .55, 
0.5). Fix the order -<o, we estimate the value of a/1 — (0.5)^ — (0.5)^ by 

(see Table m). 





(0.5,0.5)^ 


(0.5,0.5) 


E™1 sgn(fe. 


(0.5,0.5)^ 


(0.5,0.5) 



Actually - (0.5)^ - (0.5)^ = 0.7071067812. 

Example 4. Given the values of exp(x^ + ?/) at the points (2, 2.95), (2, 3.05), 
(1.95,3), (2.05,3), (1.975,2.975), (1.975,3.025), (2.025,2.975), (2.025,3.025). 

Fix the order -<o, we still choose — ^- — as our estimation 

S-l^S"(*M(2,3))-*M(2,3) 

value (see Table [5]). 

Actually exp(22 + 3) = 1096.633158. 
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i 


Yi 


exp(xf + Vi) 


interpolating value of Wj 


1 


(2, 2.95) 


1043.149728 


1043.099728 


2 


(2, 3.05) 


1152.858743 


1044.131824 


3 


(1.95, 3) 


900.0947180 


1043.425504 


4 


(2.05, 3) 


1342.783531 


1102.658424 


5 


(1.975, 2.975) 


968.3804142 


1097.459656 


6 


(1.975, 3.025) 


1018.030340 


1096.945601 


7 


(2.025, 2.975) 


1182.782509 


1096.552830 


8 


(2.025, 3.025) 


1243.425065 


1096.660126 



Table 5: estimation value of exp(2'^ + 3) 



5. Conclusion 

In this paper, we apply the Fitzpatrick algorithm to osculatory rational 
interpolation, and get the parametric solution of all the interpolation func- 
tions with the given complexity . 

For Cauchy interpolation, we present a Neville-like algorithm to determine 
the value of interpolating function at a single point without computing the 
rational interpolation function (several points can be treated similarly). Since 
each of the vectors W{i, L + 1) = {ai\y^^,hi\y^ gives an approximate value, 
we choose 



T.T=i sgn(6i 










Yo 



as our estimation value. From the examples we can see that the Neville-like 
algorithm is effective. 
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